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ABSTRACT 


The  MUSTANG,  NFS  middle  ultraviolet  spectrograph,  instrument 
was  tested  using  standard  techniques  to  determine  the  wavelength 
calibration  and  overall  sensitivity.  The  instrument  was  launched  on  March 
10,  1994  on  a  NASA  sounding  rocket  from  Poker  Flats,  Alaska.  Post-flight 
calibration  indicates  the  wavelength  calibration  did  not  change  as  a  result  of 
the  launch  and  no  significant  change  in  the  sensitivity  calibration. 
Ultraviolet  dayglow  spectra  of  the  earth's  Ionosphere  from  1800  A  to  3400 
A  were  obtained  during  a  similar  launch  on  March  19,  1992  from  White 
Sands  Missile  Range,  New  Mexico.  Data  were  obtained  on  the  downleg  of 
this  earlier  experiment  and  range  in  altitude  from  115  km  to  320  km. 
Analysis  of  the  data  from  2420  A  to  2490  A  was  conducted  to  obtain  the 
intensity  profile  of  the  011  2470.4  A  multiplet  The  analysis  used  synthetic 
spectra  generated  for  the  N2  Vegard-Kaplan  and  the  nitric-oxide  gamma 
band  emissions. 
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Z.  INTRODUCTION 


The  variaJaility  of  the  ionospheric  electron  density 
needs  to  be  better  understood  because  it  strongly  influences 
high  frequency  (HF)  radio  communications,  over-the-horizon 
(OTH)  radar,  ballistic  missile  early  warning  systems,  and 
the  Ground  Wave  Emergency  Network  (GWEN) .  In  addition  to 
this,  the  Joint  Chiefs  of  Staff  have  prioritized  the 
measurement  of  the  ionospheric  electron  density  as  number 
five  of  50  critical  global  parameters  for  inclusion  in  a 
defense  environmental  satellite  (MJCS  154-86  dated  March  21, 
1986) . 

The  ionosphere  has  the  ability  to  reflect  transmitted 
radio  waves.  The  altitude  at  which  these  radio  waves 
reflect  is  dependent  on  the  electron  density  of  the 
ionosphere  and  the  wavelength  of  the  transmitted  radio  wave. 
This  atmospheric  'mirror'  is  not  conqpletely  smooth,  but  has 
many  imperfections.  Because  of  these  imperfections  auid  their 
effect  on  radio  communication  and  radar  it  is  important  to 
be  able  to  measure  the  electron  density  profile  and  how  this 
profile  varies  over  the  entire  globe.  Currently,  measurement 
of  the  ionospheric  electron  density  is  limited  to  20  ground- 
based  ionosonde  stations  scattered  throughout  the  world. 
Although  this  is  accurate  locally,  this  method  is  not 
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capable  of  being  used  on  a  global  basis.  What  is  needed  is  a 
low-power  orbiting  system  capable  of  passively  measuring  the 
electron  density. 

Joint  NPS/NRL  research  is  currently  testing  the  Middle 
Ultraviolet  Spectrograph  (MUSTANG)  and  the  High  Resolution 
Airglow  and  Aurora  Spectrograph  (HIRAAS)  for  this  purpose. 

By  measuring  the  intensity  of  atmospheric  emissions,  the  ion 
densities  may  be  determined.  The  electron  density  may  then 
be  inferred  from  the  ion  density.  MUSTANG  was  first  launched 
aboard  a  NASA  sounding  rocket  in  March  of  1990  from  the 
White  Sands  Missile  Range,  New  Mexico,  and  again  in  a 
similar  launch  in  March  of  1992.  Analysis  of  the  data 
obtained  from  the  dotmleg  of  the  March  1992  launch  is  the 
topic  of  this  thesis. 

A.  TBSSIS  OBJECTIVES 

This  thesis  has  two  main  goals;  calibration  of  the 
MUSTANG  instrument  for  the  March,  1994  laimch  from  Poker 
Flats,  Alaska  auid  the  determination  of  the  intensity  profile 
for  the  011  2470  A  multiplet  from  data  obtained  from  the 
March  19,  1992  launch  from  White  Seunds  Missile  Range,  New 
Mexico.  Accurate  calibration  of  the  instrument  prior  to 
lavinch  is  critical  for  analyzing  the  data  obtained  during 
the  rocket  flight.  The  ultraviolet  dayglow  spectra  used  to 
determine  the  intensity  profile  were  obtained  on  the  down 
leg  of  the  1992  flight  and  range  from  115  km  to  320  km. 
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B.  raiSIS  OUTLINB 

This  thesis  is  divided  into  six  Chapters  and  one 
Appendix.  A  brief  introduction  of  the  ionosphere  along  with 
the  objectives  of  this  thesis  are  given  in  Chapter  I. 

Chapter  II  discusses  the  ionosphere  in  greater  detail  and 
briefly  introduces  atomic  and  molecular  emissions.  Chapter 
III  describes  the  instrument  in  detail  and  the  procedure  for 
collecting  the  data.  Chapter  IV  develops  the  method  and 
theory  of  the  MUSTANG  calibration  and  contains  the  final 
calibration  results.  Chapter  V  describes  the  use  of 
synthetic  spectra  and  the  method  used  in  analyzing  the  Oil 
2470  A  emission.  Finally,  Chapter  VI  concludes  with  a 
summary  of  the  results. 


ZZ.  BACKGROUND 


A.  THB  ZONOSPHERB 

The  depth  of  penetration  of  solar  radiation  incident  on 
the  earth's  atmosphere  is  dependent  on  wavelength.  Figure  1 
shows  the  altitude  of  unit  optical  depth  for  radiation  with 
wavelengths  up  to  3500  A  and  the  atmospheric  constituent 
mainly  responsible  for  its  absorption.  Also  noted  are  the 
ionization  limits  for  N2#  N,  NO,  O,  and  O2-  Unit  optical 
depth  is  the  altitude  where  the  incoming  radiation  intensity 
is  decreased  by  1/e  and  is  in^ortant  because  this  is  where 
most  of  the  energy  is  deposited.  Radiation  with  wavelengths 
longer  than  approximately  3100  A  penetrates  to  the  surface 
of  the  earth.  Radiation  from  approximately  1340  A  to  3100  A 
does  not  reach  the  surface  of  the  earth,  but  instead 
dissociates  O2  and  ozone  (O3) .  Radiation  with  wavelengths 
shorter  than  approximately  1340  A  penetrates  to  different 
altitudes,  depending  on  wavelength,  before  ionizing  the 
different  atmospheric  constituents.  Notice  the  ionization 
limits  for  all  the  major  atmospheric  constituents  are 
included  in  this  range  of  wavelengths.  This  ionizing  of  the 
atmosphere  is  the  creation  of  the  ionosphere. 

The  ionosphere  begins  cd^out  50  km  and  extends  up  to 
about  1000  km.  It  is  primarily  created  by  atmospheric 
constituents  being  ionized  by  incoming  solar  radiation  with 
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wavelengths  less  than  1340  A.  The  ionization  rate  in  the 
ionosphere  is  proportional  to  the  photon  flux  multiplied  by 
the  number  density.  Figure  2  shows  schematically  how  the 
electron  production  rate  follows  the  photon  flux  and 
particle  number  density  curves.  At  high  altitudes  there  is 
an  abiindance  of  photons,  but  few  particles.  The  reverse  is 
true  for  low  altitudes,  there  is  an  abundance  of  particles, 
but  few  photons;  most  having  been  absorbed  at  higher 
altitudes.  This  resulting  ionization  curve  is  referred  to  as 
Cha^meui  Profile.  Figure  3  is  a  more  realistic  view  of  the 
ionosphere  electron  density  profile  along  with  the 
atmospheric  constituents  contributing  to  that  density. 


Figure  1:  Unit  optical  depth  for  each  wavelength. 
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Figure  2:  Electron  production  curve. 


0+ 


Electron  Number  Density  (elect/cm''3  ) 

Figure  3;  Ionosphere  electron  density  profile.  The  solid 
line  is  electron  density  and  the  dashed  lines  are  ion 
densities . 
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The  ionosphere  is  comprised  of  four  layers:  D,  E,  Fi  and 
F2 •  The  electron  density  and  altitude  of  each  layer  depends 
primarily  on  the  intensity  of  incoming  solar  radiation.  The 
D-region  extends  from  50  km  to  90  km.  This  is  the  lowest 
region  of  the  ionosphere  and  is  created  mainly  by  short 
wavelength  ultraviolet  and  x-rays  with  the  maximum  electron 
density  occurring  edxjut  75  to  80  km.  This  layer  is  very 
weakly  ionized  with  elecron  densities  on  the  order  of  10^  - 
10^  electrons  per  cubic  centimeter.  At  night,  with  no 
incident  radiation,  the  electrons  and  ons  recombine  and  the 
D  layer  disappears.  In  D,  E,  and  FI  layer  the  primary  loss 
mechanism  is  dissociative  recombination,  for  example  NO'*'  + 
e“  -►  N  +  O.  The  major  ion  in  the  D  layer  is  NO'*'. 

The  E  layer  lies  between  90  and  140  km  with  higher 
electron  densities  ranging  from  10^  to  10^  elec/cm^ .  The 
electron  density  in  this  layer  weakens  at  night,  but  this 
layer  is  still  capable  of  reflecting  radio  waves  back  to 
earth.  The  major  ions  in  this  layer  are  N0+  and  02'*'. 

The  FI  layer  lies  between  140  and  200  km  with  electron 
densities  slightly  higher  than  the  E  layer.  The  major  ions 
in  this  layer  are  NO'*',  02'’',  and  to  a  lesser  extent  O'*'. 

The  F2  layer  extends  from  200  km  and  beyond.  In  this 
region  O'*"  is  the  major  positive  ion.  The  electron  density 
profile  peaks  in  this  region  with  densities  on  the  order  of 
2x10®  elec/cm^ ,  This  peak  is  formed  from  O'*"  diffusing 
downward  to  approximately  250  km  where  it  is  chemically  lost 
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and  the  density  begins  to  decrease.  The  primary  loss 
mechanism  is  O'*'  +  m  O  +  m'*’,  where  m  is  NO  or  O2.  The 
electron  density  in  this  layer  decreases  by  about  a  factor 
of  three  at  night,  but  this  layer  does  not  disappear.  The 
atmospheric  density  is  so  low  at  these  altitudes  that  it 
requires  more  than  one  night  for  the  ions  and  electrons  to 
recombine.  The  downward  diffusion  of  0'*'  also  plays  a  role  in 
maintaining  this  layer  at  night. 

As  mentioned  in  the  Introduction,  the  ability  to  measure 
these  electron  densities  is  important  because  of  the 
profound  effect  they  have  on  communications  and  radars. 

B.  ATOMIC  AMD  MOLBCDLAR  EMISSIONS 

Atoms  and  molecules  can  exist  in  many  discrete  energy 
states  and  often  jump  from  one  energy  state  to  another. 
Although  a  transition  from  a  higher  energy  state  to  a  lower 
energy  state  may  occur  spontaneously,  a  transition  from 
lower  energy  to  higher  energy  may  not.  Some  mechanism 
capable  of  transferring  energy  to  the  atom  or  molecule  is 
required.  There  are  many  mechanisms  capable  of  transferring 
energy,  but  this  study  is  primarily  concerned  with  only  two; 

1.  Absorption  of  an  incident  photon,  or  photoexcitation,- 

2 .  Collisional  excitation  with  an  incident  photo¬ 

electron. 

Absorption  of  incoming  solar  radiation  will  excite 
atmospheric  atoms  and  molecules  to  higher  energy  states.  If 
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this  incident  radiation  is  sufficiently  energetic,  it  will 
ionize  the  atom  or  molecule  and  create  a  photo-electron. 
These  photo-electrons  can  excite  other  atoms  and  molecules 
through  collisional  excitation.  These  excited  atoms  and 
molecules  remain  in  higher  energy  states  for  a  short  period 
of  time  before  returning  to  the  ground  state.  When  this 
transition  occurs  a  photon  with  an  energy  equal  to  the 
change  in  energy  is  emitted.  The  important  point  being  that 
each  species  of  atoms  or  molecules  can  only  transition 
between  discrete  energy  states.  Measuring  the  energy  of  the 
emitted  photon  Ceui  then  be  used  to  determine  which  species 
it  originated  from.  Since  wavelength  is  much  easier  to 


measure  than  energy,  the  relationship  ^  provides  the 
link  required.  The  MUSTANG  instrument  was  designed  to 
measure  intensity  versus  wavelength. 

In  the  case  of  the  atom,  these  transitions  primarily 
involve  a  chcuige  of  state  of  the  outer  orbiting  electrons. 
For  molecules,  the  situation  is  slightly  different. 
Molecules  can  also  rotate  and  vibrate  about  a  central  axis. 
In  addition  to  electronic  transitions,  a  molecule  will  also 
have  vibrational  and  rotational  transitions.  A  typical 
energy  level  diagram  is  shown  in  Figure  4.  As  seen  in  the 
figure,  for  any  given  electronic  transition  there  are  many 
possible  vibrational  transitions,  and  for  any  vibration 
transition  there  are  many  possible  rotational  transitions. 
The  typical  energies  for  these  transitions  are: 
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Energy 


1.  Electronic,  10  eV; 

2.  Vibrational,  .1  eV; 

3.  Rotational,  .001  eV, 


Figure  4:  Molecular  energy  level  transition  diagram. 
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III.  THE  EXPERIMENT 

A.  INTRODUCTION 

Previous  work  suggests  that  ion  densities  may  be 
determined  through  measurement  of  emissions  from  neutral 
constituents  (Bosserman,  1989  and  Danczyk,  1989) .  The 
ionospheric  electron  density  can  then  be  inferred  from  the 
ion  densities.  To  test  this  theory  NFS'  MUSTANG  and  NRL's 
HIRASS  were  launched  on  a  NASA  sounding  rocket  in  March  1990 
and  again  in  March  1992.  During  each  flight  a  door  on  the 
payload  section  opened  at  approximately  100  km  exposing  the 
two  instruments  to  the  atmosphere.  After  outgassing  the 
attitude  control  system  maneuvered  the  payload  to  the 
specified  instrument  observation- zenith-angle  (OZA) .  For  the 
1992  la\mch  the  OZA  for  the  upleg  was  100  degrees.  Data  were 
taJcen  from  150  km  to  320  km.  At  apogee  the  payload  was  again 
maneuvered  so  the  OZA  was  90  degrees.  Data  were  taken  on  the 
downleg  from  320  km  to  115  km. 

B.  INSTRUMENT  DESCRIPTION 

The  MUSTANG  instrument  is  a  l/8th  meter  Ebert-Fastie 
spectrograph  with  a  l/8th  meter  telescope  and  a  wavelength 
coverage  from  1800  A  to  3400  A.  This  wavelength  coverage  is 
the  middle  ultraviolet  (MUV)  and  part  of  the  near 
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ultraviolet  (NUV)  portion  of  the  electromagnetic  spectrum.  A 
schematic  diagram  of  the  instrument  is  shown  in  Figure  5. 


Figure  5;  Mustang  schematic  diagram.  The  major  components 
shown  are:  1)  l/8th  m  telescope  mirror,  2)  Spectrograph 
entrance  slit,  3)  l/8th  m  Ebert  mirror,  4)  Diffraction 
grating,  5)  ITT  image  intensifier,  6)  Hamamatsu  image 
detector . 


Light  which  enters  the  telescope  is  collimated  by  a 
series  of  baffles  prior  to  reaching  the  telescope  mirror. 
This  light  is  then  collected  by  a  l/8th  meter  spherical 
mirror  and  focused  onto  a  5mm  by  140  urn  vertical  slit.  After 
passing  through  the  entrance  slit  the  light  strikes  the 
Ebert  mirror.  This  mirror  collimates  the  light  and  directs 
it  onto  a  reflective  diffraction  grating  with  a  ruling 
density  of  1200  lines/mm.  The  grating  reflects  the 
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collimated  light  back  to  the  Ebert  mirror  where  it  is  then 
focused  onto  the  image  intensifier  at  the  spectrograph  focal 
plane. 

The  image  intensifier  consists  of  a  quartz  input  window, 
a  cesium  telluride  (CsTe)  photo -cathode,  two  micro  channel 
plates  (MCP)  in  cascade,  an  aluminum  screen  coated  with 
phosphor,  and  a  fiber  optic  window.  The  purpose  of  the  image 
intensifier  is  to  convert  the  ultraviolet  photons  to  visible 
photons  with  enough  gain  so  they  can  be  detected  by  the 
image  sensor.  This  is  accon^lished  when  ultraviolet  photons 
at  the  exit  focal  plane  of  the  spectrograph  strike  the 
photo-cathode  generating  photo-electrons.  These  electrons 
are  then  accelerated  toward  the  micro  chcuuiel  plates  by  an 
accelerating  voltage  of  approximately  200  Volts.  The  voltage 
across  the  micro  chamnel  plates  is  2000  Volts.  The  photo¬ 
electrons  strike  the  walls  of  the  MCP  causing  an  electron 
avalanche.  This  results  in  approximately  40,000  electrons 
being  produced  for  every  one  emitted  by  the  photo -cathode. 
These  electrons  are  again  accelerated  by  6000  Volts  between 
the  MCP  output  and  the  phosphor  screen.  The  screen 
fluoresces  as  a  result  of  the  incoming  electrons.  This 
visible  light  is  transmitted  through  a  fiber  optic  window 
with  the  same  spatial  configuration  as  the  original  spectrum 
incident  on  the  photo -cathode. 

The  photodiodes  of  the  Hamamatsu  S2300-512F  Plasma- 
Coupled  Device  (PCD)  Linear  Image  Sensor  are  sensitive  to 
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light  in  the  4000  A  to  lOOOO  A  region.  Photons  emitted  from 
the  image  intensifier  phosphor  screen  are  sensed  by  the  PCD 
image  sensor.  The  output  window  of  the  image  intensifier  and 
the  PCD  image  sensor  window  are  both  made  of  a  fiber  optic 
material  and  are  placed  in  physical  contact. 

The  light  sensitive  section  of  the  image  sensor  is  a 
25mm  by  5mm  photosensitive  area  consisting  of  512 
photodiodes  arranged  in  a  linear  array.  The  photodiodes  are 
36um  wide,  5mm  tall,  with  a  spacing  of  50um  between  the 
centers  of  each  photodiode.  The  photodiodes  convert  the 
optical  energy  of  incident  photons  to  a  charge  which  is 
proportional  to  the  corresponding  light  intensity.  The  shift 
register  transfers  this  accumulated  charge  from  the 
photodiodes  to  the  output  with  one  conqplete  spectrum  of  512 
voltages  in  series  transmitted  every  51.2  ms.  This  charge 
then  undergoes  a  10 -bit  analog- to-digital  conversion  and  is 
displayed  on  the  con^uter  or  tremsmitted  via  rocket 
telemetry. 
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IV.  CALIBRATION 


A.  NAVBLXNGTH  CALIBRATION 

The  wavelength  response  was  the  first  instrument 
calibration  to  be  performed.  This  calibration  is  required  to 
precisely  determine  what  wavelength  each  pixel  is  sensing. 
The  wavelength  calibration  was  accon^lished  using  two 
sources  with  well  known  emissions.  The  first  source  was  a 
platinum  hollow  cathode  lamp.  This  source  was  chosen  because 
it  had  23  useful  spectral  emissions  in  the  1800  A  to  3400  A 
wavelength  range.  This  lanqp  was  used  in  conjunction  with  a 
Harrison  6522A  DC  power  supply.  The  lan^  was  placed  on  the 
optical  table  with  the  power  supply  voltage  set  at  400  Volts 
amd  the  current  at  0  mA.  The  current  was  then  gradually 
increased  to  15  mA.  The  average  of  25  spectra  were  taken 
using  'Lcdsview'  and  the  groiind  support  equipment.  This 
average  is  shown  in  Figure  6.  The  23  emissions  were 
individually  isolated  and  the  centroid  of  each  pixel 
position  was  estimated. 

After  MUSTANG  was  returned  from  integration  testing,  a 
'hot  spot*  was  observed  on  the  detector.  With  no  input  to 
the  instrument  and  the  high  voltage  on,  an  abnormally  high 
output  appeared  as  shown  in  Figure  7.  The  source  of  this  is 
still  unknown.  This  odsnormal  output  was  coming  from  pixels 
134-136  as  shown  in  Figure  8.  Because  of  this,  the  platinum 
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gure  7:  'Hot  spot*  with  high  voltage  on. 
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Figure  8 :  Enlarged  view  of  ' hot  spot ' 


The  same  procedure  was  followed  for  the  second  source, 
which  was  a  mercury  lamp.  This  source  was  chosen  because  it 
added  several  emissions  in  the  longer  wavelength  range  and 
one  in  the  shorter  wavelength  of  the  detector.  The  results 
of  this  calibration  are  shown  in  Table  I. 


TABLE  I.  EMISSIONS  USED  FOR  WAVELENGTH  CALIBRATION 


PQIMT 

SOURCE 

LJLil 

PIXEL 

1 

Pt 

1916.08 

21.0 

2 

Pt 

2128.61 

87.5 

3 

Pt 

2144.23 

92.3 

4 

Pt 

2174.67 

101.6 

5 

Pt 

2202.22 

111.0 

6 

Pt 

2245.52 

125.0 

7 

Pt 

2289,27 

137.0 

8 

Pt 

2357.10 

159.6 

9 

Pt 

2440.06 

185.5 

10 

Pt 

2467.44 

195.0 

11 

Pt 

2487.17 

201.4 

12 

Pt 

2515.58 

210.0 

13 

Pt 

2539.20 

217.6 

14 

Pt 

2628.03 

246.0 

15 

Pt 

2702.40 

270.0 

16 

Pt 

2733.96 

280.0 

17 

Pt 

2830.30 

310.6 

18 

Pt 

2877.52 

325.4 

19 

Pt 

2893.86 

331.3 

20 

Pt 

2929.79 

343.2 

21 

Pt 

2997.97 

365.0 

22 

Pt 

3042.64 

379.4 

23 

Pt 

3064.71 

386.7 

1 

Hg 

1942.27 

29.0 

2 

Hg 

2536.52 

216.4 

3 

Hg 

2653.69 

254.0 

4 

Hg 

2967.28 

355.5 

5 

Hg 

3341.48 

477.0 

The  wavelengths  were  then  plotted  versus  the  pixel 
position  8uid  using  the  IDL  £\inction  Poly_£it,  a  second 
degree  polynomial  was  fit  to  the  data  with  the  following 
equation: 

X  =  -lil3xl0-^ A’* +3.1996A’  + 1850.17  A  (4-1) 

where  X  is  the  pixel  number.  This  is  shown  in  Figure  9. 
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This  equation  may  now  be  used  to  determine  the  wavelength 
associated  with  each  pixel. 


Finally,  to  test  the  accuracy  of  the  calculated  equation 
the  difference  between  the  calculated  emission  value  and  the 
actual  emission  value  was  calculated  using  the  pixel  values 
in  Table  I.  The  result  is  shown  in  Figure  10.  The  difference 
between  actual  euid  calculated  appears  to  be  random  with  only 
two  points  varying  more  than  2  A.  With  the  wavelength 
difference  between  each  pixel  being  3.133  A,  the  calculated 
equation  is  quite  accurate. 

A  second  degree  polynomial  was  used  to  fit  the  line 
because  it  produced  the  best  fit.  The  difference  between  a 
linear  fit  and  the  actual  data  is  shown  in  Figure  11  and  has 
fluctuations  about  the  zero  point  up  to  -8  A.  Clearly,  the 
second  degree  polynomial  yields  a  better  fit.  A  third  and 
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fourth  degree  polynomial  were  also  fit  to  the  data.  The 
difference  between  the  third  degree  fit  and  the  actual  data 
is  shown  in  Figure  12.  As  can  be  seen  in  the  figure,  there 
was  no  significant  improvement  for  a  third  degree  or  greater 
polynomial  fit. 


o  i oo  aoo  300  Aoo  soo 
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Figure  10;  Difference  between  actual  and  calculated  for  a 
second  degree  fit. 


Figure  11:  Difference  between  actual  and  calculated  for  a 
linear  fit. 
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Figure  12 :  Difference  between  actual  and  calculated  for  a 
third  degree  fit. 


B.  DATA  COLLECTION 

The  leJaoratory  Ground  Support  Equipment  (GSE)  provides 
the  same  timing  signals  that  the  rocket  telemetry  delivers 
to  the  instrument  during  flight  and  also  collects  the  data 
from  the  MUSTANG  instrument.  The  GSE  is  controlled  by  a 
Macintosh  computer  equipped  with  National  Instruments  Direct 
Memory  Access  (DMA)  ,  Multi -fvinction  Input /Output  (MIO)  and 
Digital  Input/Output  (DIO)  esq^ansion  cards.  Timing  signals 
and  data  acquisition  are  performed  by  'Labview',  also  by 
National  Instruments.  Power  is  provided  to  the  instrument  in 
the  laboratory  by  a  Hewlett-Packard  6206  DC  Power  Supply. 
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The  GSE  operating  procedure  is  detailed  below: 

1.  Ensure  the  two  ribbon  cables  are  connected  from  the 
Macintosh  to  GSE  box,  one  37-pin  cable  from  GSE  box 
to  MUSTANG,  two  banana  cables  from  power  supply  to 
GSE  box,  and  the  9-pin  and  15-pin  cables  from  the  5V 
and  15V  converter. 

2.  Turn  on  Macintosh  and  find  the  'MUSTANG  GSE'  file. 
Open  this  file  and  double  click  on  the  'MAIN  GSE' 
icon  to  launch  Labview  and  the  MUSTANG  program. 

3.  Click  the  r\in  button  to  start  the  program. 

4.  Click  on  the  initialize  I/O  Boards  icon.  This 
initializes  the  data  acquisition  boards  in  the 
Macintosh  computer. 

5.  Click  on  the  Initialize  Telemetry  Clocks  icon.  This 
ensures  that  the  data  acquisition  boards  provide  the 
correct  telemetry  signals  to  the  MUSTANG  instrument. 

6.  Turn  on  the  power  supply  (28  Volts) . 

7.  Turn  on  the  amp  switch  and  high  voltage  switch  on  the 
GSE  box.  The  amp  switch  turns  on  the  linear  array 
detector  and  the  high  voltage  switch  turns  on  the 
image  intensifier.  This  energizes  the  MUSTANG 
instrument . 

8.  Click  on  the  Data  Acquisition  icon.  This  launches  the 
data  acquisition  program.  It  is  important  that  step  7 
be  executed  before  step  8.  If  the  MUSTANG  instrument 
is  not  powered  it  will  not  return  a  signal.  The  data 
acquisition  program  will  wait  for  signals  in  an 
infinite  loop. 

9.  Click  on  the  Reacquire  to  get  a  new  spectrum. 

10.  Click  on  the  save  icon.  This  saves  the  spectrum.  A 
new  spectrum  can  be  generated  by  clicking  on  the 
Reacquire  icon. 

11.  Click  on  the  Quit  icon.  This  stops  the  data 
acquisition  subroutine  and  returns  to  the  GSE  panel 
window.  Now  click  on  the  Quit  icon  on  the  main  GSE 
Panel  Window.  This  stops  the  Labview  program. 

12.  Turn  off  the  power  supply  and  turn  off  the  amp  and 
high  voltage  switches  on  the  GSE  box. 

C.  SENSZTIVZTy  CALIBRATION 

1 .  Overview 

The  output  of  the  instrument  is  512  analogue 
voltages  which  are  expressed  in  decimal  units  per  pixel . 
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In  order  to  convert  this  output  in  units  of  relative 
intensity  to  an  absolute  intensity  (Rayleighs/Angstrom)  for 
use  in  measuring  the  atmospheric  emission  rate,  the  raw  data 
must  be  multiplied  by  a  conversion  factor.  The  sensitivity 
calibration  was  performed  to  determine  this  factor  which  is 
the  product  of  the  screen  reflectance  and  the  irradiance 
from  a  known  light  source  divided  by  the  instrument  output 
obtained  from  these  sources.  This  will  be  shown  below.  This 
factor  is  actually  a  function  of  wavelength  and  is  the  so- 
called  ' calibration  curve ' . 

To  obtain  this  curve  two  sources  of  known  irradicuice 
were  used.  The  EG&G  Gamma  Scientific  deuterium  lamp,  serial 
number  0155,  was  used  to  cover  the  wavelength  range  of  2000 
A  to  2700  A  and  an  Optronics  Laboratories,  Inc.,  tungsten 
filament  quartz  halogen  lamp  (FEL-IR) ,  serial  number  F-309, 
was  used  to  cover  from  2700  A  to  3400  A  range.  The 
calibrations  of  both  lan^s  are  traceable  to  the  National 
Institute  of  Science  and  Technology  (NIST) .  Because  the 
calibration  was  performed  in  open  air,  wavelengths  below 
2000  A  are  affected  by  absorption  due  to  oxygen.  Therefore, 
no  sensitivity  was  calculated  in  this  region. 

The  first  requirement  was  to  convert  the  irradiance 
of  the  sources  to  a  more  useful  unit.  The  irradiance  data 
provided  for  both  lamps  from  the  manufacturer  was  in  terms 
of  spectral  irradiance  [uW/cm2*nm]  at  50  cm.  Because  the 
output  of  the  detector  is  determined  by  the  number  of 
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photons  that  strike  each  photodiode,  the  output  of  the  two 
calibration  sources  had  to  be  converted  to  a  photon  flux. 
Using  a  spreadsheet  the  manufacturer's  irradiance  values 
were  converted  to  a  photon  flux  by  dividing  by  the  photon 
energy  at  each  wavelength.  This  resulted  in  irradiance 
values  in  units  of  (ph/cm2*s*A) .  These  irradiance  values 
were  then  multiplied  by  the  ratio  (50/distance) ^  because  the 
manufacturer  data  for  each  source  were  calculated  at  50  cm. 
These  points  were  then  fit  with  a  curve.  It  required  two 
equations  to  fit  the  deuterium  data  at  200  cm  using 
Poly_fit,  they  are: 

2400)  =  433953181?  -189935391+4.7551685x10*  (4-2) 

^  2400)  =  -0.0495457841^  +  546.1 14341^  - 2082569.01+23060629x10’  (4 -3 ) 

in  (ph/cm2*s*A]  at  200  cm  and  1  in  Angstroms. 

The  FEL-IR  source  is  considered  a  'black-body' 
radiator  with  a  temperature  of  approximately  3000  K.  To 
obtain  the  best  fit  to  the  data  we  used  a  ten^erature  of 
T=3008  K  in  the  following  eqation  for  a  black-body  radiator: 

^  2nch ^  1 

(4-4) 

where  c»3xl0® (m/s) ,  h»6 .626xl0"34 ( joule*sec) ,  k=l . 381x10*23 
(joule/k) .  This  equation  yields  irradiance  values  in 
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[uW/cin2*nm]  using  A,  in  meters.  Then,  converting  this  to 
(ph/cm2*s*A)  at  300  cm: 

Inch  1  10"’ *0.028 •A 

(4-5: 

yields  the  equation  used  to  fit  the  FEL-IR  irradiance  data 
at  300  cm.  These  curves  are  shown  in  Figures  13  and  14  for 
the  D2  and  FEL-IR  lamps,  respectively. 


Figure  13:  Curves  fit  to  Deuterium  lamp  data  at  200  cm. 


25 


These  calculations  produced  irradiance  values  for 
the  sources  in  units  of  (photons/  cm2*s*A)  at  the  respective 
distances.  These  values  are  tabulated  in  Tables  II  and  III. 


TABLE  II.  IRRADIANCE  FROM  DEUTERIUM  SOURCE 


Wavelength (Al 

Irradiance 

Energy (xlO- 

Irradiauice 

(uW/cni2*nni) 

19, 

(ph/cm2*s*A) 

at  50  cm 

( J/ph) 

at  200  cm 

2000 

4.246x10-2 

9.9320 

2.6933x10® 

2100 

4.027 

9.4620 

2.6813 

2200 

3.83 

9.0310 

2.6718 

2300 

3.69 

8.6430 

2.6897 

2400 

3.542 

8.2825 

2.694 

2500 

3.006 

7.9512 

2.3818 

2600 

2.566 

7.6400 

2.1159 

2700 

2.02 

7.3550 

1.8861 

2800 

1.908 

7.0960 

1.694 

2900 

1.647 

6.8510 

1.5145 

3000 

1.431 

6.6210 

1.3616 

3100 

1.253 

6.4080 

1.2319 

3200 

1.087 

6.2080 

1.1031 
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3300 

9.426x10-3 

6.0200 

9.8644x107 

3400 

8.457 

5.8430 

9.1184 

3500 

7.499 

5.6760 

8.3234 

TABLE  III.  IRRAIDANCE  FROM 

FEL-IR  SOURCE 

Wavelength (A) 

Irradiance 
(uW/cm2*nm) 
at  50  cm 

Energy (xlO- 
19) 

(J/ph) 

Irradiance 

{ph/cm2*s*A) 

at  300  cm 

2500 

1.87x10-2 

7.9512 

6.5853x10’^ 

2600 

3.26 

7.6400 

1.1947x10® 

2700 

5.41 

7.3550 

2.0595 

2800 

8.49 

7.0960 

3.35 

2900 

1.28x10-1 

6.8510 

5.2313 

3000 

1.85 

6.6210 

7.8235 

3100 

2.60 

6.4080 

1.136x10® 

3200 

3.56 

6.2080 

1.6056 

3300 

4.76 

6.0200 

2.2139 

3400 

6.22 

5.8430 

2.9806 

3500 

7.98 

5.6760 

3.9365 

For  MUSTANG  to  record  the  intensity  spectrum  for 
each  lamp,  the  lan^s  were  placed  one  at  a  time  in  a  closed 
box  with  an  adjusted>le  aperature.  The  lamp,  reflecting 
screen  and  MUSTANG  were  all  raised  to  the  same  height  as  the 
aperture.  The  instrument  line-of -sight  was  45  degrees  from 
the  normal  to  the  reflecting  screen.  The  normal  to  the 
screen  was  in  line  with  the  laxap  optical  axis.  The  led)  was 
then  darkened,  the  reflecting  screen  was  fully  illuminated 
by  adjusting  the  box  aperature  auid  the  spectra  were  taken 
using  the  GSE.  This  geometry  is  shown  in  Figure  15.  The 
instrument  output  for  each  lamp  is  shown  in  Figures  16  and 
17  for  D2  and  FEL_IR,  respectively. 
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Figure  16:  MUSTANG  output  for  Deuterium  lamp. 
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Figure  17;  MUSTANG  output  for  FEL-IR  lamp. 


2.  Calibration  Gaonatry  and  Thaory 

Several  characteristics  of  the  calibration  were 
established  previously  (Anderson, 1990;  Chase, 1992)  and  will 
be  assumed  in  this  calibration.  First,  it  is  assumed  that 
the  response  of  the  MUSTANG  instrument  changes  linearly  with 
changes  in  irradiauice.  This  linearity  of  response  is 
accurate  up  to  approximately  940  decimal  units/pixel 
(Chase, 1992) .  Second,  it  is  assumed  that  with  a  viewing 
angle  of  45*,  the  reflecting  screen  is  Lambertian  and  no 
correction  to  the  data  is  necessary.  Finally,  the  field  of 
view  (FOV)  of  the  instrument  was  established  as  0.06* 
horizontal  FOV,  auid  2.0*  vertical  FOV. 

The  output  of  the  MUSTANG  instrument  is  proportional 
to  the  number  of  photons  entering  the  instrument.  If  the 
output  and  the  number  of  photons  entering  is  known,  the 
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constant  of  proportionality,  'F,  can  be  determined. 

TcU^ulated  above  is  the  number  of  photons  leaving  the  source. 
Now  what  is  needed  is  the  number  of  photons  from  the 
reflecting  screen  entering  the  instrument.  This  is  done  by 
dividing  the  reflecting  screen  into  elemental  units  of  area 
and  then  integrating  over  the  area  viewed  by  the  instrument. 
Following  this,  the  objective  is  to  determine  this  photon 
flux  in  terms  of  instrument  constants  and  the  tabulated 
source  irradiance.  The  instrument  output  is  then  equated  to 
this  photon  flux  and  solved  for  in  terms  of  the  source 
irradiance,  screen  reflectivity  and  the  proportionality 
constant.  Since  the  instrument  output  is  known,  the  constant 
can  then  be  solved  for. 

The  calibration  parameter,  T,  was  defined  and 
derived  previously  (Anderson,  1990;  Chase,  1992) .  The 
following  definitions  will  be  used  in  this  discussion: 
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-spectral  irradiance 
-spectral  irradiance 


striking  the  screen; 


f  ^  ] 

Vcm*  *s*aJ 

the 


screen; 

Px  -photon  flux 


received  at  the  entremce  slit  over 


a  one -A  wide  wavelength  bin; 

At  -area  (cm^)  of  the  telescope  aperature; 

As  -area  (cm^)  of  the  spectrograph  entrance  slit; 

Ay  -area  (cm2)  of  the  screen  viewed  by  the  instrument; 
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P  -focal  length  (cm)  of  the  telescope  collecting 
mirror ; 

d  -distemce  (cm)  from  the  Ebert  mirror  to  the  screen; 

Pi  -reflectance  of  the  screen; 

I  -atmospheric  emission  rate  (R/A) ; 

r  -instrument  gain  (decimal  xinits/photon* second)  ; 
-bandwidth  of  each  pixel  (A/pixel) ; 

D  -instrument  output  (decimal  \init s/pixel)  . 

For  simplicity  the  X-siibscripts  will  be  deleted; 
spectral  irradiance/radiauice  euid  the  photon  flux  will  be 
implied  unless  otherwise  stated. 

The  geometry  of  the  calibration  setup  is  shown  above 
in  Figure  15.  Geometric  optics  allows  for  the  calculation  of 
the  area  of  the  screen  viewed  by  the  instrument,  A^,  in  the 
limit  where  the  slit  is  narrow  and  This  allows  the  use 

of  plane  wave  assumptions.  It  can  be  shown  that  the  solid 
angle  viewed  by  the  telescope  (i.e.  the  instrument  field-of 
view)  is:  fl  =  (v4, /p^).  Similarly,  the  area  of  the  screen 
viewed  by  the  slit  is  given  by:  Q=(A^cos9  /d^) .  Equating 
these  two  solid  angles  and  solving  for  gives: 


This  yields  the  area  of  the  screen  viewed  by  the  instrument. 
Next,  the  irradiance  leaving  the  reflecting  screen  is 
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needed.  For  a  Lambertian  surface,  an  ideal  diffuse 
reflector,  the  radiance  at  an  angle  0  from  the  normal  of 
the  screen  is  equal  to  the  product  of  the  irradiance  at 
0=0*  times  the  cosine  of  0  : 

l(0)  =  i.  COS0.  (4-7) 

The  area  of  the  screen  viewed  by  the  instrument  must 
be  divided  into  elemental  units  of  area  (da) .  The  solid 
angle  subtended  by  the  telescope  aperture  at  each  element  of 
area  (da)  on  the  screen  is:  The  number  of 

photons  per  second  per  euigstrom  received  at  the  telescope 
mirror  due  to  an  element  of  area  is :  dP-  L(Q  yoda . 

Substituting  for  <d  yields:  dP^*  L(dXA,  td^)da .  The  total 
number  of  photons  received  at  the  entrance  slit  of  the 
spectrograph  from  the  reflecting  screen  may  be  calculated  by 
integrating  over  the  area  of  the  screen  that  is  viewed  by 
the  instrument.  This  is  just  the  projected  area  of  the  slit 
on  the  screen.  Ignoring  reflective  losses  from  the  mirror, 
and  assuming  that  L(0  )  is  constant  over  the  elemental  area 
in  question: 


The  integral  is  equal  to  Av 


(4-8) 
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Hence : 


P=L{Q)-^A, 


(4-9) 


Substituting  Equation  (4-6)  into  (4-9)  yields; 


P=i(9)^| 


/I  ^ 


3  ^  *0)80 


A,. 


(4-10) 


Rearranging  the  terms  euid  simplifying  leaves: 


/»  = 


(4-11) 


For  a  Lambertian  surface,  where  X(9)=:Ir, cos0  : 


(4-12) 


Thus,  the  number  of  photons  per  second  per  Angstrom 
passing  through  the  entrance  slit  of  the  spectrograph  is  a 
product  of  constants  (for  constant  radiance)  euid  is 
independe>it  of  angle.  For  this  reason,  the  angle  of  the  line 
of  sight  of  the  instrument  from  the  normal  to  the  reflecting 
screen  is  not  important. 

The  spectral  radicuice  of  the  screen  is  determined  by 
equating  the  spectral  irradiance  striking  the  screen  to  the 
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total  spectral  radiance  leaving  the  screen  divided  by  the 
reflectivity  of  the  screen: 


£  = 


Ji(e)d!Ci 

i 

P 


(4-13) 


where  d!Qssin0dB<^  ,  0  ranges  from  0  to  n/2,  and  ^  ranges  from 
0  to  271.  This  forms  a  hemisphere  centered  on  the  center  of 
the  reflecting  screen.  Assuming  that  the  reflectance  is 
independent  of  6,  Equation  (4-13)  becomes: 


£  =  ^  cos8  ) 


(4-14) 


The  solution  of  this  integral  is:  EawLo/p.  Therefore,  the 
radiance  reflected  normal  to  the  screen  is: 


4 


(4-15) 


The  directional  spectral  radiance  can  now  be  rewritten  as: 


p£ 

Z(0 )  = - COS0 

n 


(4-16) 


Substituting  Equation  (4-15)  into  Equation  (4-12)  gives: 


(4-17) 
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This  is  the  photon  flux  per  Angstrom  passing  through  the 
entrance  slit. 

Converting  the  raw  data  (decimal  units/pixel)  to  a 
standard  unit  (R/A)  is  done  by  equating  the  fltix  entering 
the  spectrograph,  P,  to  the  instrument  output,  D  as  follows: 

D=prBy,;  (4-18) 

where  F  is  the  instrument  gain,  a  product  of  the  various 
efficiencies  and  transfer  functions  of  the  individual 
components  of  the  instrument,  and  By,  is  the  bandwidth  of 
each  pixel.  See  the  Appendix  for  a  definition  and  derivation 
of  r.  Substituting  for  the  flux.  Equation  (4-17),  yields: 

D  =  Tb\^^PE.  (4-19) 

The  constants  in  Equation  (4-19)  are  gathered  together  to 
define  a  new  constant: 

which  allows  the  output  of  the  instruraent  to  be  written  in 
terms  of  the  source  irradiance  and  constants: 


Now,  assuming  that  the  atmospheric  emission  rate  is 
omnidirectional  or  47il  (R/A) ,  the  spectral  radiance  received 
by  the  MUSTANG  instrument  is: 


f  10‘/) 

"I  471  J' 


(4-22) 


(  V*-  1 


where  L  is  in  units  of  — —  .  The  spectral  radiance 

\cm  sAstrJ 

recieved  by  MUSTANG  may  also  be  written  in  terms  of  the 
spectral  irradiance  of  the  atmosphere.  For  normal  incidence 
this  is: 


L-E/tc.  (4-23) 

With  omnidirectional  emissions  the  instrument  is  at  normal 
incidence  to  the  atmospheric  airglow  regardless  of  the  OZA. 
This  removes  any  angular  dependence  from  the  calculations. 
Therefore,  the  output  of  the  instrument.  Equation  (4-19) , 
due  to  the  spectral  radiance  received  from  the  atmosphere 
is: 


(4-24) 
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This  is  rewritten  in  terms  of  the  omnidirectional  emission 


rate.  Equation  (4-22) : 


(4-25) 


Solving  for  the  emission  rate  of  the  atmosphere  yields: 


/  = 


1  Y 


D; 


(4-26) 


where  I  is  in  units  of  R/A.  This  atmospheric  emission  rate 
is  now  in  terms  of  instrument  constants  and  the  output.  So, 
given  the  output  of  the  instrument,  the  atmospheric  emission 
rate  can  now  be  determined. 

A  new  constant  is  defined: 


Wbja,aJ 


(4-27) 


This  is  simplified  by  using  the  definition  for  K  from 
Equation  (4-20),  which  yields: 


Sxabstituting  this  constant  into  Equation  (4-25)  ,  where  E  is 
the  spectral  irradiance  from  a  calibration  source,  yields: 

The  instrument  output  is  now  defined  in  terms  of  the  source 
irradiance,  the  reflectivity  of  the  screen  and  the 
proportionality  constant.  The  constant,  'P,  may  then  be 
defined  in  terms  of  the  screen  reflectance,  spectral 
irradiance  and  the  response  of  the  instrument  as: 


(4-30) 


in  iinits  of  (Rayleighs/Angstrom) / (decimal  units/pixel).  This 
results  in  the  calibration  parameter,  'P,  for  the  MUSTANG 
instrument . 

Finally,  the  spectral  emission  rate,  which  will  be 
used  below,  for  the  atmosphere  is  equal  to  the  calibration 
parameter  multiplied  by  the  detector  response: 

I»D'P.  (4-31) 

This  constant,  'P,  is  dependent  on  wavelength  and  is  the 
calibration  curve  needed  to  convert  from  decimal  units/pixel 
to  Rayleighs /Angstrom.  With  the  irradiance  values  calculated 
above  and  the  instrument  output  already  obtained,  the  only 
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necessary  input  remaining  is  the  reflectance,  p,  of  the 
screen . 

3.  Reflectivity  of  screen  as  a  function  of  wavelength 

The  first  step  in  calculating  the  screen  reflectance 
curve  is  to  take  the  ratio  of  the  instrument  output  from  the 
reflecting  screen  to  the  instrument  output  looking  directly 
at  the  lamp.  This  was  done  for  each  lamp  individually. 
Because  of  the  sensitivity  of  the  instrument,  looking 
directly  at  each  lamp  would  con^letely  saturate  the  output. 
To  obtain  useable  output,  the  instrument  entrance  was 
covered  with  a  piece  of  cardboard  with  a  tiny  pinhole 
aperature.  This  reduced  the  incoming  photon  flxix  to  useable 
levels . 

Next,  the  relative  reflectance  curve  for  each  lamp 
was  normalized.  This  was  accon^lished  using  information  from 
the  mauiufacturer  concerning  the  reflectance  of  the  Barium 
Sulfate  surface  coating  of  the  screen.  At  wavelengths 
longward  of  3200  A,  this  reflectsuice  is  considered  constcuit 
at  0.97.  The  FEL-IR  ratio  curve  was  then  normalized  to  0.97 
by  multiplying  by  0.366.  The  Deuterium  ratio  was  then 
multiplied  by  0.70764  to  acheive  the  same  normalization. 
These  two  individual  curves  were  then  overlapped  to  arrive 
at  the  reflectance  curve.  The  useable  portion  of  the 
deuterium  curve  was  from  pixel  47  to  350  (2000  A  to  2951  A) . 
The  usable  portion  of  the  FEL-IR  curve  was  from  pixel  250  to 


450  (2640  A  to  3259  A) ,  this  yielded  an  overlap  of  311  A. 
Outside  of  these  ranges  the  normalized  curves  were  far  too 
noisy.  This  left  the  ranges  of  1850  to  2000  A  and  3259  to 
3400  A  unaccounted  for.  For  the  3259  to  3400  A  range,  the 
reflectance  was  assumed  constant  at  0.97.  For  the  1850  to 
2000  A  range  a  line  was  fit  to  the  deuterium  lamp  data  with 
the  equation: 


y  =.001983x+034  (  4  -  32) 

where  x  is  the  pixel  number.  This  final  reflectance  curve  is 
shown  in  Figure  18. 

With  the  reflectance  curve  now  determined,  the 
calibration  parameter,  'F,  could  now  be  calculated  using: 

,4.33) 

This  calibration  curve  was  calculated  individually  for  each 
lamp  as  shown  in  Figures  19  and  20.  The  two  curves  were  then 
overlayed  and  merged  at  pixel  206,  2499.73  A,  to  yield  the 
final  calibration  curve  as  shown  in  Figure  21. 
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Figure  18:  Reflectance  curve. 
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Figure  19 : Deuterium  lamp  sensitivity  curve. 
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D.  POST-PLZOBT  CALZBSATION 

Post -flight  calibration  indicated  that  the  wavelength 
calibration  did  not  change  and  there  was  no  significant 
change  in  the  sensitivity  calibration  as  a  result  of  the 
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launch.  Figure  22  shows  a  comparison  of  the  pre-flight  and 
post-flight  platinum  spectra.  For  wavelength  calibration, 
the  position  of  each  peak  with  respect  to  pixel  position  is 
the  important  factor,  not  intensity.  As  can  be  seen,  the 
pixel  position  of  each  emission  peak  is  essentially 
unchanged,  indicating  no  change  in  the  wavelength 
calibration.  Figure  23  is  the  post-flight  deuterium  spectrum 
compared  with  the  pre-flight  deuterium  spectrum.  Again,  the 
pre-flight  spectrum  is  plotted  directly  over  the  top  of  the 
post -flight  spectrum  indicating  no  chauige  in  the  sensitivity 
calibration. 


Figure  22:  Post-flight  platinum  spectrum  plotted  over  the 
pre-flight  platinum  spectrum. 
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post -flight  deuterium  spectrum. 
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V.  DATA  AKALYSZS 


A.  INTRODUCTION 

The  objective  of  this  analysis  is  to  estimate  the 
temperature  profile  from  the  NO  y(0,1)  and  Vegard-Kaplan 
(0,5)  bands,  and  the  intensity  profile  for  the  singly 
ionized  oxygen  (Oil)  2470  A  multiplet  on  the  downleg  of  the 
March  1992  rocket  experiment.  The  spectra  for  the  downleg 
were  taken  from  320  km  down  to  115  km.  Shown  in  Figure  24 
is  a  representative  spectrum  taken  at  180  km.  It  is  a  plot 
of  intensity  versus  wavelength. 
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Figure  24:  Representative  spectrum  taken  at  180  km. 


The  determination  of  the  2470  A  intensity  profile  is 
complicated  by  the  presence  of  contaminating  emissions  from 
the  nitric-oxide  gamma  (y)  transition  and  the  nitrogen  (N2) 
Vegard-Kaplan  (VK)  transition.  The  nitric  oxide  y  bands  are 
the  result  of  the  electronic  transition  from  the  first 
excited  state  to  the  ground  state  (A^L^-^X^ri)  and  extend 
from  1873  A  to  6126  A.  As  stated  ed}ove,  within  this  one 
electronic  transition  are  mcuiy  possible  vibrational 
transitions.  For  example,  y(0,6}  refers  to  the  transition 
from  the  zero  vibrational  level  of  the  first  excited  state 
to  the  sixth  vibrational  level  of  the  ground  state.  The  NO  y 
bands  are  produced  primarily  by  photo-excitation. 

The  molecular  nitrogen  (N2)  Vegard-Kaplan  bands  are  also 
the  result  of  the  transition  from  the  first  excited  state  to 
the  groxmd  state  (A^Z+^-^xlSg"*')  .  The  emissions  in  this 
transition  range  from  1250  A  to  5325  A.  There  is  a 
difference  between  this  transition  emd  the  NO  y  trouisition 
cUxjve.  This  transition  is  a  forbidden  treuisition  auid  is 
produced  by  collisions  with  photo-electrons  vice  photo¬ 
excitation. 

The  011  2470  A  emission  is  produced  by  the  transition 
from  the  second  excited  state  to  the  ground  state  (^P->^S) 
of  singly  ionized  atomic  oxygen.  Unlike  the  molecular 
transitions  described  cdjove,  the  2470  A  transition  involves 
no  vibrational  or  rotational  energy  levels  and  does  not 
produce  a  band  of  emissions.  Like  the  N2  VK  transition,  this 
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transition  is  also  forbidden  and  is  produced  by  collisions 
with  photo-electrons  vice  photo-excitation. 

B.  QBNBRATZON  OP  SniTHETIC  SPBCTRA 

The  models  generating  the  synthetic  spectra  for  N2  VK 
and  NO  y  bands  below  were  developed  by  Cleary  (1986) .  These 
synthetic  spectra  are  generated  based  on  the  quantum 
selection  xxiles  and  the  probadiilities  governing  the 
different  treuisitions .  The  probedsilities  of  the  various 
vibrational  transitions  within  each  electronic  state  are 
given  by  the  Franck-Condon  factors.  These  synthetic  spectra 
are  then  convolved  with  the  instrument  slit  functions  to 
match  the  output  of  the  instrument.  There  are  five  different 
slit  functions  for  the  MUSTANG  instrument.  The  slit  function 
that  produced  the  best  fit  was  used  in  the  convolution. 
Because  the  011  2470  A  emission  doesn't  depend  on 
temperature,  the  synthetic  spectrum  is  generated  by  simply 
inserting  a  delta  function  at  the  required  wavelength  and 
convolving  it  with  a  slit  function. 

C.  FITTZMG  THE  DATA  WITH  THE  SYNTHETIC  SPECTRA 

The  data  used  to  analyze  the  2470  A  multiplet  consisted 
of  42  records  in  5  km  increments  that  ranged  from  115  km  to 
320  km.  To  fit  each  spectrum,  synthetic  spectra  were 
generated  for  NO  y  and  N2  VK  bands.  These  synthetic  spectra 
were  multiplied  by  constants  to  match  the  actual  intensities 
recorded  by  MUSTANG.  Once  the  contribution  from  the 


contaminating  emissions  were  estimated,  they  could  be 
subtracted  out .  The  remaining  feature  is  attributed  to  the 
011  2470  A  emission.  This  feature  is  fit  by  multiplying  the 
convolved  delta  function  by  a  constant  in  the  same  manner  as 
above . 

The  spectrum  i_t  each  altitude  was  fit  individually.  The 
synthetic  models  require  an  estimated  rotational  temperature 
for  the  specified  altitude  and  the  vibrational  transition  of 
interest  for  input.  By  running  the  models  repeatedly  using 
different  '«^'»brational  transitions,  theoretical  spectra  for 
the  NO  Y  and  N2  VK  bands  were  developed.  The  NO  y  bands 
contributing  to  the  2100  A  to  2800  A  wavelength  range  are: 
(0,0),  (0,1),  (0,2),  (1,0),  (1,1),  and  (1,3).  The  N2  VK 

bands  contributing  to  the  2100  A  to  2800  A  wavelength  range 
are:  (0,3),  (0,4),  (0,5),  (0,6),  (1,4),  (2,5),  and  (3,6). 

The  on  2470  A  emission  is  contaminated  by  the  VK(3,6)  band 
and  the  y(0,2)  band.  The  synthetically  generated  spectra  for 
the  beuids  listed  abo've  are  shown  plotted  over  the  180  )cm 
spectrum  in  Figure  25. 
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bands  of  interest  plotted  over  the  180  km  spectrum. 

The  rotational  temperature  also  affects  the  shape  of  the 
synthetic  spectra.  At  higher  rotational  temperatures  more 
molecules  are  in  higher  rotational  states  causing  the  shape 
of  the  molecular  band  to  shift  as  shown  in  Figure  26.  This 
temperature  dependence  affects  NO  y  and  N2  VK  differently 
since  the  N2  VK  bands  are  degraded  to  the  red  and  the  NO  y 
bands  are  degraded  to  the  blue.  This  means  that  for  the 
synthetic  N2  VK  bands,  using  a  higher  rotational  temperature 
shifts  the  shape  of  the  molecular  band  to  longer,  or  'red' 
wavelengths.  For  the  synthetic  NO  y  bands  a  higher 
rotational  temperature  shifts  the  shape  of  the  molecular 
beuid  to  shorter  or  'blue'  wavelengths.  Since  atomic  oxygen 
emissions  involve  no  vibrational  or  rotational  energy 
states,  changing  the  temperature  has  no  effect  on  the  shape 
of  the  feature.  The  actual  atmospheric  temperature  at  the 


49 


altitude  of  interest  can  be  inferred  from  the  rotational 
temperature  that  best  fits  the  data.  This  was  performed  on 
all  42  spectra.  The  results  are  listed  in  Table  IV  and  the 
temperature  profile  is  shown  in  Figure  27.  The  rotational 
temperature  profile  can  be  modeled  using: 

Tialtitude)  =  7.  -  (7.  -  7;,,)  *  ^  ^  5  _ . 


where  s=0. 0255/km,  Tqo=1300  K  and  Tii5=550  K.  The 

temperature  at  115  km  is  considered}ly  higher  than  predicted 
by  MSIS,  the  standard  atmospheric  model.  More  work  is 
required  to  resolve  this.  A  comparison  of  these  two  profiles 
is  shovm.  in  Figure  28. 


spectrum  by  using  a  higher  rotational  temperature. 
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Figure  27;  The  estimated  temperature  profile  for  the 
downleg . 
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Figure  28;  Comparison  of  temperature  profile  inferred  from 
the  data  with  profile  calculated  from  Equation  (5-1) .  The 
solid  line  is  from  Equation  (5-1) . 


TABLE  IV. 

ESTIMATED  TEMPERATURES 

115 

550 

220 

1250 

120 

575 

225 

1250 

125 

670 

230 

1300 

130 

730 

235 

1300 

135 

810 

240 

1300 

140 

920 

245 

1300 

145 

950 

250 

1300 

150 

950 

255 

1300 

155 

1000 

260 

1300 

160 

1000 

265 

1300 

165 

1000 

270 

1300 

170 

1000 

275 

1300 

175 

1050 

280 

1300 

180 

1100 

285 

1300 

185 

1100 

290 

1300 

190 

1150 

295 

1300 

195 

1200 

300 

1300 

200 

1200 

305 

1300 

205 

1200 

310 

1300 

210 

1200 

315 

1300 

215 

1250 

320 

1300 

To  fit  the  NO  and  N2  data,  vibrational  bands  had  to  be 
chosen  that  showed  very  little  contamination  from  other  non- 
related  emissions.  By  multiplying  the  theoretical 
intensities  by  constants  to  match  these  selected  emissions 
the  contribution  from  these  molecules  to  the  2470  A  feature 
can  be  estimated  and  subtracted  out.  This  is  shown  in  Figure 
29.  The  remaining  intensity  to  this  feature  is  attributed  to 
the  011  2470  A  emission.  A  representative  final  fit  is  shown 
in  Figure  30. 
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Ficfure  29:  Estimating  contribution  to  2470  A  feature  by  NO  y 
and  N2  VK  emissions.  The  solid  line  is  the  observed 
intensity,  the  dashed  line  is  the  sum  of  NO  y  and  N2  VK. 
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Figure  30:  Representative  final  fit  for  180  km  showing 
contributions  from  all  emissions.  Curves  presented  are:  1) 
NO  y,  2)  N2  VK,  3)  011  2470  A,  4)  Overall  fit,  5)  Data. 


Because  the  nitric  oxide  y  transition  is  an  allowed 
transition  and  is  produced  by  photo-excitation,  the  volume 


emission  rates  of  these  bands  depend  on  the  intensity  of  the 
incoming  solar  radiation.  This  incoming  solar  radiation  can 
be  assumed  to  be  constant  over  the  rocket  launch  time  frame. 
The  NO  scale  factor  can  be  used  to  estimate  the  NO  column 
density  by  starting  with: 

Data(R  /  A)  =  Synthetil — —  . — 1  * 

\pfxel*  molecule*  sec)  (5-2) 

Solving  for  the  scale  factor  yields: 


Data 

Synthetic 


=  Scalt 


R/A 


ph!  pix*  mol*  sec 


} 


10‘ 


ph 


crn^  *  sec  t 


3.133-^ 

pvc 


(5-3) 


which  reduces  to: 


Data 

Synthetic 


3.133x10* 


(5-4) 


The  scale  factor  multiplied  by  a  constant  yields 
column  density.  These  values  for  the  downleg  are 
in  Table  V  and  plotted  versus  altitude  in  Figure 


the  NO 

tabulated 

31. 
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Figure  31:  Estimated  nitric  oxide  column  density. 


TABLE  V.  ESTIMATED  NITRIC  OXIDE  COLUMN  DENSITY 


Alt  (Ian) 

Column  Density 

115 

(em~2) 

3.76el5 

120 

2.66 

125 

2.06 

130 

1.72 

135 

1.44 

140 

1.22 

145 

1.06 

150 

1.00 

155 

8.61el4 

Alt  (km) 

Column  Density 

(cm~2) 

220 

1.50 

225 

1.47 

230 

1.34 

235 

1.25 

240 

1.09 

245 

1.03 

250 

8.61el3 

255 

8.45 

260 

7.83 

Because  the  N2  VK  and  the  011  2470  A  are  forbidden 
transitions,  the  volume  emission  rates  are  dependent  on  the 
photo- electron  flux.  This  flux  is  not  constant  and  varies 
greatly,  so  the  scale  factors  for  the  N2  VK  band  and  the 
2470  A  feature  are  used  to  estimate  the  intensities  vice  the 
column  densities.  The  total  intensity  is  solved  for  by 
integrating  to  find  the  area  under  the  molecular  band  or  the 
area  under  the  2470  A  feature  as  shown  in  Equation  (5-5) : 

(5-5) 


But,  this  intensity  can  be  estimated  by  using  the  scale 
factor  and  the  synthetically  generated  spectral  emission  as 
shown  in  Equation  (5-6) : 


/ »  \syr(Xi* Scale* dk^ Scale* &X*Y,Syn . 


(5-6) 


where  13 3  (A/pixel)  and  is  the  total  intensity  or 

area  under  the  synthetically  generated  spectral  emission 
after  it  has  been  convolved  with  the  slit  function  and 
before  it  is  multiplied  by  the  scale  factor.  The  total 
intensity  of  the  synthetically  generated  N2  VK(0,5)  band 
with  a  temperature  of  550  K  is  7.4xl0"2.  The  difference 
between  the  intensity  at  550  K  and  1300  K  is  only  a  few 
percent.  As  a  result,  this  total  intensity  was  used  for  all 
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altitudes.  Because  the  synthetically  generated  spectrum  for 
the  2470  A  feature  is  a  delta  function,  the  area  is  1.0. 


This  calculation  for  the  intensity  was  performed  for 
each  spectrum.  The  resulting  intensity  values  for  the  N2 
VK(0,5)  band  are  listed  in  Table  VI,  with  the  intensity 
profile  shown  in  Figure  32.  The  resulting  intensity  values 
for  the  2470  A  feature  are  listed  in  Ted>le  VII,  with  the 
intensity  profile  shown  in  Figure  33. 
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:  011  2470  A  intensity  profile. 
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TABLE  VI. 

ESTIMATED  NITROGEN 

(N5)  intensities 

AIL  (tan) 

Intenaity  (R1 

Ait  (tan) 

Intensity  (n) 

115 

2090.64 

220 

1509.91 

120 

2253.24 

225 

1300.84 

125 

2137.10 

230 

1254.38 

130 

2439.08 

235 

1091.78 

135 

2555.22 

240 

905.94 

140 

2787.52 

245 

836.25 

145 

3019.81 

250 

708.49 

150 

3484.40 

255 

638.80 

155 

3252.10 

260 

557.50 

160 

3600.54 

265 

441.35 

165 

3670.23 

270 

452.97 

170 

3368.25 

275 

360.05 

175 

3437.94 

280 

325.21 

180 

2926.89 

285 

308.95 

185 

2833.98 

290 

220. SI 

190 

2555.22 

295 

197.44 

195 

2555.22 

300 

164.92 

200 

2322.93 

305 

134.73 

205 

2206.78 

310 

109.17 

210 

1881.57 

315 

116.47 

215 

1649.28 

320 

69.68 
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TABLE  VII. 

ESTIMATED  Oil 

2470 

A  INTENSITIES 

JEL 

115 

1567.97 

220 

2289.24 

120 

1662.05 

225 

2257.84 

125 

1567.97 

230 

2289.24 

130 

1348.46 

235 

2257.88 

135 

1505.25 

240 

2257.88 

140 

1473.89 

245 

2257.88 

145 

1662.05 

250 

2195.16 

150 

1787.49 

255 

2195.16 

155 

1787.49 

260 

2132.44 

160 

1693.41 

265 

2069.72 

165 

1724.77 

270 

1881.57 

170 

1850.21 

275 

1912.93 

175 

2038.36 

280 

1850.21 

180 

1975.65 

285 

1787.49 

185 

2101.08 

290 

1818.85 

190 

2038.36 

295 

1630.69 

195 

2007.00 

300 

1630.69 

200 

2195.16 

305 

1536.61 

205 

2195.16 

310 

1348.46 

210 

2289.24 

315 

1379.82 

215 

2289.24 

320 

1003.50 

In  order  to  verify  the  aOaove  results,  the  intensity  of 
the  N2  VK(0,5)  band  was  estimated  directly  from  the  data 
using  an  integration  approximation  technique.  The  shape  of 
the  VK(0,5)  molecular  band  can  be  approximated  by  a 
triangle.  Using  area* (1/2) *base*height,  the  intensity  can  be 
estimated.  This  is  shown  schematically  in  Figure  34.  Using 
this  method,  the  area  under  the  N2  VK(0,5)  band  at  180  km 
was  estimated  to  be  2400  (R) .  Whereas,  using  the  fitting 
method  edx>ve  the  intensity  was  estimated  to  be  2600  (R) ,  the 
difference  between  these  values  being  only  7%. 
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Figure  34:  Use  of  estimation  for  total  intensity. 
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VI.  CQNCLUSZON 


A.  SUmORY  OP  FZKDINOS 

The  work  in  this  thesis  involved  two  distinct  aspects. 
First,  the  MUSTANG  instrument  was  calibrated  in  preparation 
for  the  March  1994  launch  from  Poker  Flats,  Alaska.  This 
calibration  included  a  wavelength  and  a  sensitivity 
calibration.  The  wavelength  calibration  fit  a  second  degree 
polynomial  to  the  wavelength  versus  pixel  plot  of  the 
platinum  and  mercury  lamp.  This  polynomial  is  used  to 
determine  the  wavelength  of  each  pixel  in  the  image  sensor. 
The  sensitivity  calibration  calculated  the  calibration  curve 
needed  to  convert  the  instrument  output  to  an  absolute 
intensity  in  Rayleigh's/Angstrom.  This  calibration  also  used 
two  lamps  of  known  output  and  required  determining  the 
reflectivity  of  the  reflecting  screen  prior  to  calculating 
the  final  curve.  The  instrument  was  launched  on  March  10, 
1994  and  post-flight  tests  indicated  that  no  significant 
change  occurred  in  the  calibration  as  a  result  of  the  rocket 
launch. 

The  second  area  of  research  was  the  estimation  of  the 
temperature  profile  from  the  NO  y(0,l)  and  VK(0,5)  beuids, 
and  the  determination  of  the  intensity  profile  of  the  Oil 
2470  A  multiplet.  Data  from  the  downleg  of  the  March  1992 
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rocket  launch  from  the  White  Sands  Missile  Range  in  New 
Mexico  consisted  of  42  records  ranging  in  altitude  from  115 
km  to  320  km.  By  fitting  synthetically  generated  spectra  for 
NO  y(0,l)  and  N2  VK(0,5)  to  the  actual  data,  the 
contribution  from  these  two  emissions  to  the  2470  A  feature 
could  be  estimated  and  siabtracted  out.  The  remaining 
intensity  is  attributed  to  the  011  2470  A  emission  and  was 
fit  with  a  synthetically  generated  spectrum  similar  to  the 
procedure  above.  The  constants  used  to  fit  the  synthetic 
data  to  the  actual  data  were  then  used  to  estimate  the 
nitric  oxide  column  density,  the  N2  VK(0,5)  intensity 
profile,  the  011  2470  A  intensity  profile  and  the 
temperature  profile.  These  results  will  now  be  used  as 
inputs  to  other  models  to  infer  the  electron  density  profile 
for  115  km  to  320  km. 

B.  RBCOMMBllDATZON  FOR  FURTHER  RESEARCH 

During  the  calibration  a  'hot  spot'  was  discovered. 
Further  research  into  the  cause  of  this  ' hot  spot '  and  the 
effect  it  has  on  the  output  of  the  instrument  would  be 
beneficial. 
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APPENDIX 


The  instrument  gain,  F,  in  decimal  units  per  photon  per 
second  is  the  product  of  the  various  efficiencies  and 
transfer  functions  of  the  individual  components  of  the 
instrument.  In  order  to  derive  this  quantity  the  following 
definitions  are  needed: 


P  -photon  flux 


^s*AJ 


received  at  the  entrance  slit  over 


a  one -A  wide  wavelength  bin 
111  -quantum  efficiency  (el/ph)  of  the  phot o- cathode ; 
Fmcp  -gain  (el /el)  of  the  microchannel  plates; 
r\2  -efficiency  (ph/el)  of  the  phosphor  screen; 

T]3  -quantum  efficiency  (el/ph)  of  the  photo  diode; 

T  -integration  time  of  photodiode  (sec) ; 
gi  -treuisfer  function  (V/el)  of  the  image  detector 
amplifier; 

g2  -transfer  function  (decimal  units/V)  of  the  analog- 
to-digital  conversion; 

By,  -bandwidth  of  each  pixel  (A/pixel)  . 
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Over  one  integration  time,  the  output  from  a  photodiode 
(Wp<j)  will  be: 

(A-1) 

in  electrons  per  Angstrom. 

During  this  same  period  of  time,  the  instrument  output  will 
be : 

=  (A-2) 

Substitute  Equation  (A-1)  into  Equation  (A'-2)  ,  and 
rearranging  terms  gives: 

D  *  P(glg2ni^mq,n2n3f  )K  •  (A-3  ) 

Therefore,  the  overall  MUSTANG  instrument  gain,  F,  may  be 
defined  as: 

r  *  )  •  (A-4  ) 
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